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HAMILTONIANS BY BOGOLIUBOV TRANSFORMATIONS 
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Abstract. We provide general conditions for which bosonic quadratic 
Hamiltonians on Fock spaces can be diagonalized by Bogoliubov trans¬ 
formations. Our results cover the case when quantum systems have infi¬ 
nite degrees of freedom and the associated one-body kinetic and paring 
operators are unbounded. Our sufficient conditions are optimal in the 
sense that they become necessary when the relevant one-body operators 
commute. 


1. Introduction 

We consider Hamiltonians on Fock space which are quadratic in terms of 
bosonic creation and annihilation operators. In many cases the quadratic 
Hamiltonians can be diagonalized by Bogoliubov transformations, namely 
they can be transformed to those of noninteracting particles by a special 
class of unitary operators which preserve the CCR algebra. The aim of our 
present work is to give rigorous conditions for which the diagonalization can 
be carried out for quantum systems of inhnite degrees of freedom where the 
kinetic and paring operators are unbounded. 

1.1. Quadratic Hamiltonian. Let us introduce the mathematical setting. 
Our one-body Hilbert space f) is a complex separable Hilbert space with 
inner product (.,.) which is linear in the second variable and anti-linear in 
the first. In the grand canonical ensemble the number of particles is not 
fixed and it is natural to introduce the (bosonic) Fock space 

OO N 

:= 0 (g) = C © fi © (1) fl) © • • • 

N=0 sym 

Noninteracting systems are described by the Hamiltonians of the form 

OO N 

dr(/i) := ^ /ij) = 0©/i©(/i©l©l©/i)©--- 

N=o j=i 

on the Fock space, where h > 0 is a self-adjoint operator on ff. Although 
we will work in an abstract setting, the reader may keep in mind the typical 
example that h = —A+V{x) on ff = where V is an external potential 

which serves to bind the particles. The operator dT{h) is well-defined on 
the core 

M n 

U ©©"w 

M>0 n=0 sym 
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and it can be extended to a positive self-adjoint operator on Fock space by 
Friedrichs’ extension. The spectrum of dF(/i) is nothing but the closure of 
the finite sums of elements of the spectrum of h. In particular, the spectrum 
of the particle number operator M := dF(l) is {0,1,2, 

In many physical situations the interaction between particles plays a cru¬ 
cial role and it complicates the picture dramatically. In principle, solving 
interacting systems exactly is mostly unrealistic and certain approximations 
are necessary. In the celebrated 1947 paper |lj, Bogoliubov introduced an 
approximation theory for a weakly interacting Bose gas where the many- 
body system is effectively described by a quadratic Hamiltonian on Fock 
space, which will be described below. We refer to the book m for a ped¬ 
agogical introduction to Bogoliubov’s approximation. Bogoliubov’s theory 
has been justified rigorously in various situations including the ground state 
energy of one and two-component Bose gases Ha dg [25], the Lee-Huang- 
Yang formula of dilute gases [ElliniET] and the excitation spectrum in the 
mean-field limit [sziiiiiiiaiaiSQ]. 

To describe quadratic Hamiltonians on bosonic Fock space, we introduce 
the creation and annihilation operators. For any vector / G f), the cre¬ 
ation operator a*{f) and the annihilation operator a{f) are defined by the 
following actions 

a*ifN+l) j fa(l) ® ... (8) L(N) j = ^ Y ® ••• ® 

\<t&Sn j CTeSjv+1 

a{fN+l) Y ® ••• ®/f7(Af) j = ViV Y {fN+lJa(l))fa{2)®-®fu(N), 

for all /i,...,/Ar+i in f), and all N = 0,1,2,.... These operators satisfy the 
canonical commutation relations (CCR) 

[a{f),a{g)] =0, [a*{f),a*{g)] = 0, [a{f),a*{g)] = {f,g), V/, 5 Gf). (1) 

In particular, for every / € 1) we have 

a (/)| 0)=0 

where |0) = 1 0 0 0 0 ... is the Fock space vacuum. 

In general, a quadratic Hamiltonian on Fock space is a linear operator 
which is quadratic in terms of creation and annihilation operators. For 
example, dF(/i) is a quadratic Hamiltonian because we can write 

dT{h) = ^ (/m,^/n)o*(/m)a(/n) 
m,n>l 

where {fn}n>i 0 D(h) is an arbitrary orthonormal basis for [) (the sum on 
the right side is independent of the choice of the basis). In this paper, we 
will consider a general quadratic operator of the form 

m = dT{h) + ^ Y {{J*kfm,fn)a{U)a{fn) + {J*kfmJn)a*{fm)a*{fn)) ( 2 ) 

m,n>l 
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where A: :[)—)•{)* is an unbounded linear operator with D{h) C D{k) (called 
pairing operator) and J : 1) —1)* is the anti-unitary operator 0 defined by 

= {f,g), V/,5rEf). 

Since H remains the same when k is replaced by {k + Jk*J)/2, we will always 
assume without loss of generality that 

k* = J*kJ*. (3) 

In fact, the formula m is formal but H can be defined properly as a 
quadratic form as follows. For every normalized vector 'I' E .F(f)) with finite 
particle number expectation, namely ('I',AA'h) < oo, its one-particle density 
matrices yijj : f) —f) and a-q/ : f) —f)* are linear operators defined by 

{f,l^g) = {^,o*{g)a{f)^), {J f,a^,g) = {^,a*{g)a*{f)^), yf,gGi). (4) 
A formal calculation using ([2]) leads to the expression 

+ ^TT{k*a^). (5) 

The formula ([5]) makes sense when and k*aq/ are trace class 

operators. We will use ([5|) to define El as a quadratic form with a dense 
form domain described below. 

Since f) is separable and D{h) is dense in f), we can choose finite-dimensional 
subspaces {Qn}’^=i such that 

Q1CQ2C ... C D{h) and [J = I). 

n>l 

Then it is straightforward to verify that 

M 

e^= U U (6) 

M>0n>l sym 

is a dense subspace of T'(f)). Moreover, for every normalized vector T E 
Q, 7 ii( and are finite rank operators with ranges in D{h) and JD{h), 
respectively. Note that JD{h) C D{k*) because D{h) C D{k) and k* = 
J*kJ*. Thus and k*aii, are trace class and ('I',]HI4') is well- 

defined by ([5]) for every normalized vector 'I' in Q. 

Under certain conditions (see Lemma [9]), we can show that the quadratic 
form El defined by ([5]) is bounded from below and closable, and hence its 
closure defines a self-adjoint operator by m Theorem VIII. 15]. Let us 
mention that if k is not Hilbert-Schmidt, then the vacuum does not belong 
to the operator domain of El (although it belongs to the form domain). 
Therefore, it is not easy to define H as an operator with a dense domain at 
the beginning. 

A key feature of the quadratic Hamiltonians in Bogoliubov’s theory is 
that they can be diagonalized to those of noninteracting systems by a special 
class of unitary operators which preserve the CCR algebra. By Bogoliubov’s 
argument [3], the diagonalization problem on Fock space can be associated 
to a diagonalization problem on f) 0 f)* in a very natural way which will be 
briefly recalled below. 

^If C : I) —>■ it is anti-linear, then C* : .S —>■ I) is defined by {C*g, /)(, = {Cf, g)si for all 
/ G fog G .ft. The anti-linear map C is an anti-unitary if C*C = R and CC* = 1^;. 
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1.2. Bogoliubov transformation. Since we will consider transformations 
on f) 0 1^*, it is convenient to introduce the generalized annihilation and 
creation operators 

Mf ® Jg) = a{f) + a*{g), A* {f ® Jg) = a*(f) + a{g), V/, 5 Gf). (7) 

They satisfy the conjugate and canonical commutation relations 


A*{Fi) = A{JFi), [A{F^),A*{F2)] = {Fi,SF2), VFi,T 2 G f) 0 ^ 


where we have introduced the block operators on 1} 0 f)* 


5 = 




J 


0 J* \ 

J 0 J ' 


( 8 ) 

( 9 ) 


Note that S = S~^ = S'* is a unitary on f) 0 1^* and J' = = J'* is 

an anti-unitary. The symplectic matrix S is the bosonic analogue to the 
identity in the fermionic case. 

We say that a bounded operator V on f) 0 ()* is unitarily implemented by 
a unitary operator Uv on Fock space if 


VvA{F)Vl; = A{VF), VTGf)0l}*. (10) 

It is easy to see that if (fTOl) holds true, then the OCR ([8]) imply the following 
compatibility conditions 


JVJ = V, V*SV = S = VSV*. (11) 


Any bounded operator V on f) 0 1)* satisfying m is called a Bogoliubov 
transformation. See [6l Chap. 11] for an alternative description of Bogoli¬ 
ubov transformations in the context of symplectic geometry. 

The condition JVJ = V means that V has the block form 


V 


U J*VJ* \ 
V JUJ* J 


( 12 ) 


where ?7 ; f) —)• f) and R ; f) —)• [)* are linear bounded operators. Under this 
form, the condition V*SV = S = VSV* is equivalent to 


U*U = l + V*V, UU* = 1 + J*VV*J, V*JU = U*J*V. (13) 


It is a fundamental result that a Bogoliubov transformation V of the form 
m is unitarily implementable if and only if it satisfies Shale’s condition 

m 


IHS = Tr(U*U) < oo. (14) 

Now we come back to the problem of diagonalizing H. Using the formal 
formula O) and the assumptior 


where 


A:= 


= J*kJ*, we can write 



(15) 

h k* \ 


k JhJ* ) 

(16) 


E {Fm,AFn)A*{Fm)A{Fn). 

m,n>l 


and 
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Here {Fn}n>i is an orthonormal basis for f) 0 1^* and the definition B [_4 is 
independent of the choice of the basis. Note that JAJ = A because of the 
symmetry condition k* = J*kJ*. 

Now let V be a Bogoliubov transformation on 1)01)* which is implemented 
by a unitary operator Uv on Fock space as in (11011 . Then we can verify that 
UvHyiUv ~ Hvyiv* and hence (fT5|) is equivalent to 

UvHUt; = ^ Tr(/i). (17) 

In particular, if is block diagonal, namely 




e 0 

0 J^J* 


for some operator ^ : I) —)• f), then (fT7t) reduces to 

UvHUt; = dr(0 + ^TV(e-/i). (18) 

Note that all formulas (fT^ . (fTTll and (fTHll are formal because h, ^ and ^ — /i 
may be not trace class. Nevertheless, the above heuristic argument suggests 
that the diagonalization problem on HI can be reduced to the diagonalization 
problem on A by Bogoliubov transformations. 


1.3. Diagonalization conditions. In this paper we are interested in the 
conditions on h and k such that A and El can be diagonalized rigorously. 

Let us make some historical remarks. The physical model in Bogoliubov’s 
1947 paper [1] corresponds to the case when dim f) = 2 and .4 is a 2 x 2 real 
matrix which can be diagonalized explicitly (more precisely, in his case par¬ 
ticles only come in pairs with momenta ±p and each pair can be diagonalized 
independently). In fact, when dimi) is finite, the diagonalization of A by 
symplectic matrices can be done by Williamson’s Theorem [26]. We refer to 
Hormander |13] for a complete discussion on the diagonalization problem in 
the finite dimensional case. 

In the 1950’s and 1960’s, Friedrichs |9| and Berezin |3| gave general di¬ 
agonalization results in the case dim!) = Too, assuming that h is bounded, 
k is Hilbert-Schmidt, and .4 > /i > 0 for a constant fj,. Note that the gap 
condition A> fi requires that h> fj. > 0. 

In the present paper we always assume that .4 > 0 but we do not require a 
gap. In some cases, the weaker assumption .4 > 0 might be also considered, 
but is is usually transferred back to the strict case .4 > 0 by using an 
appropriate decomposition; see Kato and Mugibayashi m for a further 
discussion. 

In many physical applications, it is important to consider unbounded 
operators. In the recent works on the excitation spectrum of interacting 
Bose gases, the diagonalization problem has been studied by Grech and 
Seiringer in mi when /i is a positive operator with compact resolvent and 
k is Hilbert-Schmidt, and then by Lewin, Nam, Serfaty and Solovej |15l 
Appendix A] when /i is a general unbounded operator satisfying /i > /r > 0. 

Very recently, in 2014, Bach and Bru [2] established for the first time 
the diagonalization problem when h is not bounded below away from zero. 
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They assumed that /i > 0, ||A:/i“^|| < 1 and kh~^ is Hilbert-Schmidt for all 
s G [0,1 + e] for some e > 0 (see conditions (A2) and (A5) in [2]). 

In the present paper, we relax not only the gap condition h > fj, > 0 
but also the Hilbert-Schmidt conditions on k {k is even allowed to be an 
unbounded operator). Our conditions are motivated by the following simple 
example where all relevant operators commute. 


Example (Commutative case). Let h and k be multiplication operators 
on I) = L^(fI,C), for some measure space 0. Then J is simply complex 
conjugation and we can identify I)* = I) for simplicity. Assume that /i > 0, 
but k is not necessarily real-valued. Then 

h k \ ^ ^ J, 

h A ) >" . 



if and only if —1 < G < 1 with G := \k\h In this case, A is diagonalized 
by the linear operator 


-G 


V ■.= 




2Vl - G2 


-G 




1 


in the sense that 




e 0 

0 C 


with ^ := /i\/1 — = s/hF^-lA > 0. 


It is straightforward to verify that V always satisfies the compatibility condi¬ 
tions m- Moreover, V is bounded (and hence a Bogoliubov transformation) 
if and only if ||G|| = ||/c/i“^|| < 1 and in this case 

||V||~(l-||G||)->/* (19) 

(which means that the ratio between ||V|| and (1 — ||G||)“^/^ is bounded from 
above and below by universal positive constants). By Shale’s condition (fTTD . 
V is unitarily implementable if and only if kh~^ is Hilbert-Schmidt and in 
this case, under the conventional form (|12l) . 

III^IIhs ~ (1 - ||G||)“^/^||G||hs- (20) 


Finally, from (1181) and the simple estimates 

— -k^h~^ > i — h = h? — k"^ — h > —k‘^h~^ 

we deduce that IH is bounded from below if and only if kh~^^‘^ is Hilbert- 
Schmidt and in this case 


inf (t(]HI) ~ — ||A:/i 


( 21 ) 


Thus, in summary, in the above commutative example we have the fol¬ 
lowing optimal conditions: 

• .4. is diagonalized by a Bogoliubov transformation V if and only if 

< 1 . 

• V is unitarily implementable if and only if kh~^ is Hilbert-Schmidt. 

• HI is bounded from below if and only if kh~^^‘^ is Hilbert-Schmidt. 
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The main message of our work is that the above necessary and sufficient 
conditions in the commutative case are indeed sufficient also in the general 
non-commutative case. Our results are formulated precisely in Theorem [1] 
and Theorem [2] in the next section. 
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grant No. 2012/07/N/ST1/03185 and the Austrian Science Fund (FWF) 
project Nr. P 27533-N27 (MN). 


2. Main results 


In this section we state our main results and explain the strategy of the 
proof. Our first main result concerns diagonalization of block operators. 


Theorem 1 (Diagonalization of bosonic block operators). 

(i) (Existence). Let : 1} —>■ f) and A: : f) —)• f}* 6e (unbounded) linear 
operators satisfying h = h* > 0, k* = J*kJ* and D{h) C D{k). Assume 
that the operator G := h~^^‘^J*kh~^^‘^ is densely defined and extends to a 
bounded operator satisfying ||G|| < 1. Then we can define the self-adjoint 
operator 

^'■={1 jhj - ) >" 

by Friedrichs’ extension. This operator can be diagonalized by a bosonic 
Bogoliubov transformation V on f) © I)* m the sense that 


VAV* 


0 J^J* ) 


for a self-adjoint operator ^ > 0 on f). Moreover, we have 

(i + mV'" 

Vi-iiGiiy • 



( 22 ) 


(ii) (Implementability). Assume further that G is Hilbert-Schmidt. Then V 
is unitarily implementable and, under the block form m, 

||P||hs < (23) 

Next, we consider the diagonalization of quadratic Hamiltonians. 

Theorem 2 (Diagonalization of quadratic Hamiltonians). We keep all as¬ 
sumptions in Theorem [T] (that ||G|| < 1 and G is Hilbert-Schmidt) and 
assume further that kh~^^‘^ is Hilbert-Schmidt. Then the quadratic Hamil¬ 
tonian H, defined as a quadratic form by is bounded from below and 
closable, and hence its closure defines a self-adjoint operator which we still 
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denote by H. Moreover, i/Uv is the unitary operator on Fock space imple¬ 
menting the Bogoliubov transformation V in TheoremUl then 

UvHUv = dr(0 + inf ct(]H). (24) 

Finally, El has a unique ground state 'I'o = Uy|0) whose one-particle density 
matrices are 7 >ifp = V*V and aifp = JU* J*V and 

inf(T(H) = + > -]^\\kh-^/^\\\^. (25) 

In particular, and k*aq,Q are trace class. 

Before explaining the proof strategy, let us give some remarks. 

• Since 

GG* < ||GG*|| = \\Gf 

we have 

kh-^k* < \\GfJhJ*. (26) 

Thus, the boundedness condition ||G|| < 1 is a strengthened version of the 
positivity ^ > 0. This follows from the following elementary fact whose 
proof is left to the reader. 

Lemma 3 (Positivity of block operators). 

(i) Let c, d : f) —>• 1) and b : i) ^ 1)* be bounded operators such that c > 0, d > 0 
and b* = J*bJ*. Then 

> 0 on 1) © 1)* i/ and only if hc~^h* < d. 

(ii) If c, d, b are unbounded operators such that c, d are self-adjoint and d > 
bc~^h* (which in particular requires that bc~^b* is densely defined), then D > 
0 as a quadratic form on the domain D{c^/'^)®D{d}-/‘^). Consequently, D can 
be extended to a non-negative self-adjoint operator by Friedrichs’ extension. 
Moreover, if bc~^b* < d, then Z? > 0. 

• The condition ||G|| < 1 can be interpreted as a non-commutative ana¬ 
logue of the bound || A:|| < 1 in the commutative case. In general, ||G|| < 1 
is weaAer than ||A:/i“^|| < 1 because 

{kh-^k*f = kh-^k*kh-^k* = Jh{h-^k* Jh-^k*kh-^J*kh-^)hJ* 

< \\h-^k*Jh-^k*kh-^J*kh-^\\Jh^J* < \\kh-^f{JhJ*f 

and the square root is operator monotone. 

• The implementability condition ||G||hs < oo can be interpreted as a non- 

commutative analogue of the condition ||A:/i“^||hs < oo in the commutative 
case. In general, these two conditions are not comparable, but if we assume 
further that < oo then ||G||hs < oo is weaker than ||/c/i“^||hs < 

oo. Indeed, if kh~^^‘^ is Hilbert-Schmidt then kh~^k* is trace class and we 
have the spectral decomposition 

J*kh-^k*J='^Xn\Un){Un\ 
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for an eigenbasis {un} for f). Therefore, 

||G||^s = J*kh-^k* Jh-^y/ J*kh-^k* J) 

= '^Xn{Un,h~^Un) = '^{Un, J* kh~^ k* Jh~^Un) 

n n 

= Ti{rkh-^k*Jh-^) = l^{{rkh-^f) < 

Here we have used the identity Tt{W*W) = Tr(TnT*) which holds true for 
every linear operator W. 

• The condition ||A:/i“^/^||hs < oo is the same as in the commutative case. 
This necessary condition was proved by Bruneau and Dereziiiski in [5] when 
k is Hilbert-Schmidt. Note that in order to ensure that El is bounded from 
below, we do not really need the conditions ||G|| < 1 and ||G||hs < oo in 
Theorem [U see Lemma [9l 

• Our estimates (f22]) . (i2H|) and ([251) are comparable to the formulas (fT9l) . 
(j20p and (I2ip in the commutative case. In particular, the implement ability 
bound (l23l) provide an upper bound on the particle number expectation of 
the quasi-free state Uy|0) because ||H||^g = (0|UvA/’Uy|0). 

• In Bach and Bru’s work [2], they proved (|24p under the assumptions 

||/c/i“^|| < 1 and ||A:/i“*|| HS < oo for all s G [0,1 + e]. 

Our Theorem [2] is an improvement of Bach and Bru’s result and our condi¬ 
tions seem to be optimal for (|24p to hold. 

• When BI is not bounded from below, the diagonalization result on 
A in Theorem [1] is of its own interest. In particular, assuming that V is 
unitarily implementable, we can consider the positive self-adjoint operator 
Uvdr(^)Uy as a renormalized version of H (shifted by an infinite constant). 
It allows us to discuss several important properties of the physical system, 
such as the (renormalized) ground state Uy|0) and the excitation spectrum, 
although El is not bounded from below. 

Sketch of the proof. The starting point of our approach is to employ a 
connection between the bosonic diagonalization problem and its fermionic 
analogue. Such kind of connection has been known for a long time; see 
Araki [T] for a heuristic discussion. To be precise, we will use the following 
diagonalization result for fermionic block operators. 

Theorem 4 (Diagonalization of fermionic block operators). Let B be a self- 
adjoint operator on f)© f}* such that JBJ = —B and such that dimKer(H) 
is either even (possibly 0) or infinite. Then there exists a unitary operator 
U on f} © f)* such that JUJ = U and 

"^"' = (o -J°J 

for some operator ^ > 0 on f). Moreover, i/Ker(H) = {0}, then ^ > 0. 

By applying Theorem H] to B = , with S given in Q, we can 

construct the Bogoliubov transformation V in Theorem [T] explicitly: 

V :=U\B\^/‘^A-^/‘^. 
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This explicit construction is similar to the one used by Simon, Chaturvedi 
and Srinivasan [23] where they offered a simple proof of Williamsons’ The¬ 
orem. The implement ability of V is proved using a detailed study of V*V = 

It should be mentioned that when H is bounded from below and V is 
unitarily implement able, the desired formula (1241) . 

UvHUv = dr(0 + inf u(H), 

does not follow immediately. In fact, it is not easy to carry out the formal 
argument in Section 11.21 and in particular we could not prove that ^ — h is 
trace class. Instead, we will use the quadratic form expression of H and a key 
ingredient in our proof is the kinetic energy estimate < oo. 

Organization of the paper. We will prove Theorems|31[Tjand[2jin Sections 
[3l m and (5] respectively. 

3. DiAGONALIZATION of block OPERATORS: FERMIONIC CASE 

In this section we prove Theorem |31 Let us start by recalling a well-known 
diagonalization result on anti-linear operators. 

Lemma 5 (Diagonalization of anti-linear operators). Let \] he a separable 
Hilbert space and let C he an anti-linear operator on f}. Assume that C = 
C* and is either compact or equal to 1. Then C has an orthonormal 
eigenbasis for 1} with non-negative eigenvalues. 

The proof of Lemma [5] is provided in the Appendix. As a consequence, 
we have 

Lemma 6. Let Sj be a separable Hilbert space such that dim.^ is even or 
infinite. Let J = J* = be an anti-unitary on Sj. Then there is a 

subspace A of S) such that 

© JA. 

Proof. Let dim.^ = 2d be even or infinite. By Lemma [5l we can find an 
orthonormal basis {un}^=i for 1) such that Jun = Un for all n > 1 (note 
that I is the only non-negative eigenvalue of J because = I). Define 

+ iu2j-i 

■ 

Then the vectors U form an orthonormal basis for 9). 

Therefore, = .^ © J A. with A. = Span{nj}^^;^. □ 

Now we are ready to give 

Proof of Theorem^ Since B is self-adjoint, we can define the spectral pro¬ 
jection 1(B > 0) and 1{B < 0). By the spectral theorem, we can decompose 

[) © 1)* = P+ © Ker(B) © P_ 

where 

P+ := 1{B > 0)(1} © r), P_ := 1{B < 0)(li © 

The condition JBJ = —B implies that P- = ffP+ and that ff leaves 
Ker(P) invariant. Since fL = J* = J~^ is an anti-unitary on Ker(P) and 
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dimKer(i?) is even or infinite, according to Lemma[ 6 l there is a subspace ^ 
of Ker(i?) such that 

Ker(i?) = ^0 J 

Thus we have 

(P+ 0 © J{P+ 0 = fi 0 r = (f) © 0) © © 0). (27) 

Let W : © ^ ^ f) 0 0 be an arbitrary unitary (which exists because 

0 ^ and f} have the same dimension). From (I27p it follows that 

U ■.= W®JWJ 

is a unitary on 1)0f)*. It is also clear from the definition ollA that JUJ = lA. 

Now we show that lABlA* is block-diagonal. Note that for every / G 1), 
we have W*{f 0 0) G 0 and hence BW*{f 0 0) G P+ by the spectral 
property of B. Therefore, 

VFSVF*(/0O) G 1)0 0, V/Gl). 

This observation allows us to define a linear operator ^ : 1) —1) by 

{^f)(B0:=WBW*{f(B0), V/G 1). 

Note that ^ > 0 because 

if, ?/) = (/ © 0, m © 0) = {W* if 0 0),BW* if 0 0)) > 0 (28) 

for all / G 1). The last inequality follows from the facts that W*{f © 0) G 
0 and that the restriction of B on © .^ is nonnegative. 

We will now show that 

"^"• = (o -J°r) 

which is equivalent to 

lABlA*(feJg) = Uf}e(-J^g), Vf,gef). (29) 

For every / G 1), we have U*{f © 0) = W*{f © 0) G P+ © .^. Therefore, 
BW{f 0 0) = BW*{f 0 0) G P+ and hence 

UBU*{f®0) = WBW*{f®0) = {^f)®0, V/Gl). (30) 

Similarly, for every (7 G 1), we have Z7*(O0 J 5 ) = Jg) G 

Therefore, BU*{0 0 Jg) = BJW*J{0 © Jg) G P- and hence 

UBU*{Q 0 Jg) = JWJBJW*J{0 0 Jg) = -JWBW*{g 0 0) 

= -J{{^g)e0) = -JCg, V 5 GI). (31) 

Here we have used JBJ = —B. The desired formula ()29p follows imme¬ 
diately from (f30P and (f3Tp . When Ker(H) = {0}, then C > 0 because the 
inequality (f28p is strict for every / / 0 . □ 
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4. DiAGONALIZATION of block OPERATORS: BOSONIC CASE 


Proof of Theorem [IJ (i) Existence. Let us consider 


B := 


It is clear that B is self-adjoint and Ker(i?) = {0} because ^ > 0. Moreover, 
JBJ = —B because JAJ = A and JSJ = —S. By applying TheoremlH 
we can find a unitary operator W on f} 0 f)* such that JIAJ = lA and 


"^"• = (o -JfJ' 

for some self-adjoint operator ^ > 0 on f). 
Using (I26p and Lemma [3l we find that 

6A < SAS < 6-^A with 


= : D 


e l-||g|l 
i + iiGir 


Therefore, 

5A^ <B^ = A^/^SASA^^^ < 5-^A^, 


and hence 

6^/'^A < \B\ < 6-^/^A 


(32) 

(33) 


because the square root is operator monotone. Consequently, the operator 
| 5 | 1 / 2 ^- 1/2 

is well-defined on D{A) and it can be extended to a bounded 
operator on f) 0 f)*. Thus 

V :=U\B\^/^A-^/^ 


is a bounded operator on f) 0 f)* and, by (1331) . 

||V|| = < j-1/4^ 


(34) 


It is straightforward to see that V is a Bogoliubov transformation. Indeed, 
because J commutes with A, B and U, it also commutes with V. Moreover, 

= U\B^I‘^U*{fAB-^U*)U\B\^l‘^U* 

= \D\^/^D-^\D\^/^ = S 


V*5V = A-^/‘^\B\^/^U* SU\B\^/‘^A-^/^ 

= {U\B\ )S{U\B\ )UA-^/^ 

= A-^/‘^U*\D\^/‘^S\D\^/‘^UA-^l‘^ 

= A-^/‘^U*DUA-^/‘^ = A-^/‘^BA-^/‘^ = S. 


Finally, V diagonalizes A because 

VMV* = U\B\U* = |D| = ^ 

(ii) Implementability. We shall show that Shale’s condition Tr(U*U) < 
oo is equivalent to ||V*V — I||hs < cc and give explicit upper bounds on 
||V*V- 1||hs and ||U||hs in terms of ||G||hs- 
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Step 1. Our starting point is the formula 

V*V = 

From the functional calculus, we have 

dt 


H = 


1 r 

TT Jo 


(35) 


(36) 


t + H^Vi 

for every self-adjoint operator H >0. Therefore, from (|35ll it follows that 
V*V - 1 = 

B‘^ A^ 


-io 

i / 


Using the resolvent identity 

1 1 


t + B^ t + A^ 
1 1 
t + A^ ~ t + B^ 


^ :(i?2_^2^ 1 


(37) 


t + B2 


t + A^ t + B"^ t + A^ 
and the expression 

B‘^ -A^ = A^/‘^SASA^/‘^ -A^ = A^I’^EA^I^ 

where 

E = E* := SAS -A = -2 


we find that 


1 


1 


= 


0 k* 
k 0 


^ EA^/^- ^ 


t + A^ t + B^' 

can be estimated by the Cauchy-Schwarz 


t + A^ t + B^ 

for all t > 0. The right-side of 
inequality 

±{XY + Y*X*) < e-^XX* + eY*Y 

with 

X := 


(38) 


^ EA^^^\B\-^, Y ■.= \B\—K^A-^^^. 


t -|- ^2 ’ ' ^ 2?2‘ 

Note that XY = Y*X* because the left-side of (j38p is self-adjoint. There¬ 
fore, from ([38|1 it follows that 

52 ^2 


± 2A~^^‘^ 


t + B^ t + A^ 




< e 


-1 


1 


= e 


-1 


t + ^2 
1 


ejAI'^b-'^jAI'^e 


1 


t + ^2 




B^ 


(t+ .B2)2 




ESA-^SE 


1 




B^ 




t + ^2—^“t + ^2 (t + S2)2-'‘ 

for all e > 0. Here we have used B’^ = AAI'^SASA}^!’^ and S~^ = S in the 
last equality. From the functional calculus, we have 

^■2 


B = 


2 r 
TT Jo 


0 {t + H^A 


Vtdt 


(40) 
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for every self-adjoint operator H > 0. Let us integrate (f39]l against the 
measure 2'K~^y/tdt on (0,oo), then use (1371) for the left side and use ()iO]l 
and (f35]l for the right side. This gives 

±2(V*V - 1) < e-^K + + eV*V (41) 


where 


K ■.= 


-I 

TT Jo 


oo 1 1 

ESA-^SE—^^ftdt. 


t + A^ 


t + A^ 


Step 2. Now we show that K is trace class. Using the equality Tr(lU*lU) = 
Tr(inU*), which holds for every operator IT, and the formula (|4np we get 




^.ESA-^SE-^ 
A^ t + A^ 


Vidt 


Vesa-^se 


1 


(t + A^y 


VESA-^Se] Vtdt 


Tr (^ESA-^SE 0 ^ 


{t + A^Y 

= Tr (VESA-^SEA-WESA-^SE 


Vidt] VESA-^SE 


By Lemma [3] and assumption kh < ||G|pJ/iJ*, we have 


A 


1-lie'll 


> ^0 : — 


h 0 


which is equivalent to (1 — ||G||)^ ^ ^ because the inverse mapping is 


0 JhJ 

- llOrllJ^ - ^ A~^ 

operator monotone. Thus 

Tr(iL) = Tr (^VESA-^SEA-WESA-^SE 

1 


> 0 , 


< 


1 - ||G|| V 


1 


< 


1-IIGII 

1 

(1-1® 


Tr ( AQ^^^ESA-^SEAn^^^ 


Tr ( Ao^^^ESAa^SEAn^^^ 


(42) 


Here we have used Tr(lU*lU) = Tr(lUlU*) again. From the explicit formulas 


E = -2 


0 k* 
k 0 


) Aq — 


h 0 
0 JhJ* 


S = 


1 0 
0 -1 


it is straightforward to compute 
1 - 1/2 


Aq"^"ESAq^SEAq 


- 1/2 ^ 4 


G*G 0 
0 JG*GJ* 


Inserting this formula into 


and using k* = J*kJ* we obtain 


Tr(iL) < 


(1-||G||)2 


Tr(G*G) < oo. 


(43) 
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Step 3. Now we bound ||V*V — 1||hs and ||t^||HS using (HT|) and (H3|) . By 
combining (HD) and the simple bound ||V|| < (cf. (|M]) i. we deduce 

that 

±2{V*V-1) <e-^K+ (44) 

for all e > 0. We have the following general fact whose proof can be found 
in the Appendix. 


Lemma 7. Let L = L* be a bounded operator on a separable Hilbert space. 
If there exists a trace class operator K > 0 such that 

±2L <e-^K + e, Ve > 0, (45) 

then L is Hilbert-Schmidt and ||L||^g < Tr(A'). 

From (j43D . (1441) and Lemma [7| we conclude that 

||VW- l||?is < Tr(i^) < (46) 

Note that (|45D implies immediately an upper bound on ||B||hs- Indeed, from 
the block form (|1‘2I) of V and the relations (|13p we hnd that 

ffj.) (47) 

where 

X:=V*V and Y := JV*JU = JU* J*V. 


Using (fT^ again we have 

Y*Y = {JU* rvy JU* rv = v*juu*j*v = x + x‘^. (48) 

Consequently, 


||VW - lll^s > 4nils > 4Tr(A) = 4||U|||s. 
Therefore, from (1461) we can conclude that V is Hilbert-Schmidt and 

2(1 + 11011 ) 1/2 


HS ^ 


■l|G|| 


2 

HS- 


( 1 - 11011 ) 5/2 

To obtain the better bound (I23p we need to analyze (I4ip more carefully. We 
will need the following technical result. 


Lemma 8. Assume that A :[)—>■ f) is a nonnegative trace class operator 
and y : p) —>■ f}* is a Hilbert-Schmidt operator satisfying 

Y* = J*YJ*, Y*Y = X + X^. 

Then there exists an orthonormal basis {un}n>i for 1) and non-negative num¬ 
bers {An}n>i such that 

Yuyi — XfiUfi^ Yu^i — JY JUji — X^Ju^ Vtz ^ 1. (49) 

Proof of Lemma\^ Note that C := J*Y is an anti-linear operator on p) sat¬ 
isfying C = C* and = y*y = A -|- A^. Thus both A and are 
trace class and they have common eigenspaces. For every such eigenspace 
VF, since C leaves W invariant, it follows from Lemma [5] that C admits an 
orthonormal eigenbasis for W with non-negative eigenvalues. Thus there ex¬ 
ists an orthonormal basis {un}n>i for f} and non-negative numbers {A„}n>i 
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such that Xun = XnUn and Cun = \/Xn + X^Un for all n. Here we have 
used = X + X‘^ for the relation between eigenvalues of X and eigenvalues 
of C. The conclusion ()49p then follows immediately. □ 

Now let {un}n>i satisfy (09]). Using ( 07 ]) we can rewrite m as 

2(2 - e) ( y + e, Ve > 0. 

Taking the expectation against Un 0 Jun we hnd that 

4(2 — e)(An 0 0 A^) < e ^(un (B Jun, Xun (B Jun)-I-2e (50) 

for all n > 1 and for all e > 0. We will show that 

^ ^ 0 7 Xvjji 0 J , Vn^l. (^4) 

Indeed, if A,i < 1/4, then (l50l) implies that 

8-\/K - 2e < 4(2 - e)< £~^{un 0 Jun, Kun 0 Jun) + 2e, 
and hence 

< e~^{Un® JUn,KUn® Jun)+4:6, Ve > 0. 

Optimizing over e > 0 leads to 00 ). On the other hand, if A^ > 1/4, then 
we can choose e = 1/2 and deduce from (l50)) that 

1 0 8Xn < 4(2 — £){\n + \/Xn + A^) < 2{Un 0 JUn, KUn 0 JUn) + 1 

which is equivalent to (1511) . Summing (]5ip over n > 1 and using the fact 
that {2~^/‘^Un 0 Jun\n>i is an orthonormal family in f) 0 1 )* we get 

Tr(u*l/) = ^A„ < Jun,Kun® Jun) < ^Tr(i7). 

n n 

Inserting the upper bound (1431) into the latter estimate we obtain ()23p . □ 


5. Diagonalization of quadratic Hamiltonians 

In this section we will prove Theorem |7J Recall that H is dehned as a 
quadratic form by ([0 

(T,eT) = +^TT{k*a^) 


where and are the one-particle density matrices of 4' defined by 0 p. 
From the OCR ([0, we see that a* = J*aJ* and 

where the generalized creation and annihilation operators A*{F), A{F) are 
dehned by ([0. Thus 


75 , 


a 




* I — 


> 0 . 


(y.\^ 1 “h 

By Lemma O it follows that 

75 , > 0 and 75 , > q:|,(1 0 J 75 , 


(53) 


( 54 ) 
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Consequently, since 

(l + J7vpJ*)-'>(l + ||7M>ll)-' 

we have 


7^(1 + ||7'l'||) ^ Q:|,av]>. 

Furthermore, if has finite particle number expectation, namely < 

oo, then is trace class and a-qt is Hilbert-Schmidt. 

Note that from (fM)l and the assumption < oo, it follows that 

if both 7 ,j, and are trace class, then ka^, = ■ h}/‘^a\, is 

trace class (because kh~^^‘^ and are Hilbert-Schmidt), and hence the 

right side of dS]) is well-defined and finite. 

Now we start proving Theorem [5J The first step is to verify that El is 
bounded from below. 


Lemma 9 (Lower bound on El). Assume that h > 0, k* = J*kJ*, kh ^k* < 
JhJ* and Tr:{kh~^k*) < oo. Then 

H > -^TT{kh-^k*) 

as a quadratie form. Moreover, if kh~^k* < 5JhJ* for some 0 < 5 < 1, 
then 

fd fd 

(1 + Vd)dr{h) + ^ Tr{kh-^k*) > M > (1 + V6)dr{h) - ^ TT{kh-^k*) 

as quadratic forms. Consequently, the quadratic form El defines a self- 
adjoint operator, still denoted by El, such that inff 7 (EI) > —^Tv{kh~^k*). 

Note that we do not require that G = h~^/‘^ J*kh~^/‘^ is Hilbert-Schmidt 
in Lemma [9j The first lower bound in Lemma [9] was proved in [5l Theorem 
5.4] under the additional condition that k is Hilbert-Schmidt. Our proof 
below is different from [5] and it does not require the boundedness of k. 


Proof of Lemma O Let 41 be a normalized vector in the domain Q defined 
in ([ 6 ]). Since 7 >i. and aq; are finite-rank operators, we can use the cyclicity 
of the trace and the Cauchy-Schwarz inequality to write 

I Tiik*aq,)\ = I Tt{aq,k*)\ = \ Tr((l + J-fq,J*)-^/^aq,h^/^h-^/^k*{l + Jjq,J*y/^)\ 

< 11(1 + Jjq,J*)-^/^aq,h^/^ns\\h-^^^k*{l + Jjq,J*)^/^HS 

.1 1/2 


= TT{h^^^a*q,{l + J-fq,J*)-^aq,h^/^) 

X TV(^(1 +J7vi>J*)^/2^/i"^r(l +J7^J 
Using a|,(l -|- J^q,J*)~^aq, < 'jqi and kh~^k* < JhJ*, we get 


11/2 


(55) 


71/2 


Tiik*aq,)\< TT{h^/^-fq,h^/^) Tr(fc/i-ir)+Tr(/i^/V/i^^^) 

< U{h^/^-fq,h^/‘^) + i Tv{kh-^k*). 


1/2 


(56) 
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Here in the last estimate we have used ■\/x{x + y) < x+y/2 for real numbers 
x,y > 0. Thus by definition ([5|), we get 

(T,1IT) = + ^Ti{k*a^,) > -]^TT{kh-^k*) 

for T E Q. Thus M > -{l/2)Ti{kh-^k*). 

Now we assume further that kh~^k* < 6JhJ* for some 0 < (i < 1. Then 
inserting this bound into ([55]), we can improve (l56]) to 

1 


Tr(ravi>)| < V6(^TT{h^/'^-f^h^/‘^) + ^ Tr(fch-^r)). 


Therefore, 


(T, MT) = + U Tr(r ovt) 

> (1 - ^) Tr(/i^/2^M>/i^/2) _ ^ Tr{kh-^k*) 


/T 

= (1 - v^)(T,dr(/i)^) - ^ Tiikh-^k*), 


and similarly. 


/T 

(T,1IT) < (l + V5)(T,dr(/i)^) + ^Tr(M-ir) 

for T E Q. In summary, 

/T /T 

(1 + VS)dT{h) + ^ Tv{kh-^k*) > e > (1 + VS)dT{h) - ^ TT{kh-^k*) 

as quadratic forms. Hence, the form domain of HI can be extended to be 
the same with that of dr(/i), which is closed because dT{h) is a self-adjoint 
operator. Since the quadratic form H is bounded from below and closable, 
its closure defines a self-adjoint operator by |2T1 Theorem VHI.15]. The 
quadratic form bound HI > —(1/2) Tr{kh~^k*) remains valid as an operator 
bound. □ 


In order to prove the diagonalization formula ()24p and the identity 
we need to show that both hdl‘^V*Vhdl‘^ and k*JU*J*V are trace class. To 
this end, we will use the following lemma which is inspired by Grillakis and 
Machedon ini eq. (42)] (see also [IS])- Recall that A and V have the block 
forms (fT2|) andfTUl respectively. 


(57) 


Lemma 10 (Diagonalization equations). //V.4V* is block diagonal, then 
j hX -Xh + k*Y - Y*k = 0, 

j JhJ*Y + Yh + kX + JXJ*k + k = 0. 

where X = V*V and Y = JV*JU = JU*J*V. 

Note that by computing the off-diagonal term of VAV* we get 

uhv* + j*vhu*r + uk*ju*j* + rvrkv* = o (ss) 

which is essentially equivalent to dSTj). However, the equations m are 

easier to analyze because they are linear (in terms of X and Y), while (15811 

is nonlinear (in terms of U and V). 
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Proof of Lemma\JU Recall that V~^ = SV*S. Since V^V* = is 

block diagonal, VASV~^ is also block diagonal. Therefore, = 

0, and hence [^5, V^^S'V] = 0. Using V~^ = SV*S again and (H7D we have 


y-^SV = SV*V = 

Therefore, 

0 = = 


1 + 2X 2Y* 

-2Y -{1 + 2JXJ* 


h -k* 
k -JhJ* 


1 + 2X 2Y* 

-2Y -{1 + 2JXJ*' 

h -k* 


(59) 


= 2 


1 + 2X 2Y* 

-2Y -{1 + 2JXJ*) )[ k -JhJ* 

hX -Xh + k*Y - Y*k hY* + Y* JhJ* + k*JXJ* + Xk* + k* 
JhJ*Y + Yh + kX + JXJ*k + k JhXJ* - JXhJ* + kY* - Yk* 


and (f57)l follows immediately. 


□ 


Now we are ready to prove 

Lemma 11. Under conditions of Theorem\^ the operators hfl'^Xhfl'^ and 
k*Y are trace class, where X = V*V and Y = JV*JU = JU*J*V. 


Proof. By Lemma [8] we can find an orthonormal basis {rin}n>i for 1) and 
non-negative numbers {An}n>i such that 

Xuyi — XfiUn, YUn — JY Ju^i — 's/ A^T XfJun, Vn ^ 1. 

Using the second equation in (f57)l . we obtain 

0 = {Jun, {JhJ*Y + Yh + kX + JXJ*k + k)un) 

— ^'s/Xji T X^iuji, hun) T 2(1 “b 2\yi){^J kuyi ), Vn ^ 1. 

By the Cauchy-Schwarz inequality 

\{JUn,kUn)\ = \{h~^^‘^k* JUn,h^^‘^Un)\ < \\h~^/‘^ k* JUn\\ ' \\h^^'^Un\\ 

we find that 

Xn{Un,hUn) < \\h~^ /k* J Unf ■ 


for all n > 1. Summing the latter estimate over n > 1 and using 


(1 + 2Xny 

4 + AXn 


< An + — < 


All + 


1 


we obtain 


Ti{h^/^Xh^/^) = Tr{X^/^hX^/^) = E 

n 

< riAII + iWllh-l/AVllnf 

' n 

Ti{kh~^k*) < 00 . 

Thus h^/^Xh^/^ is a trace class operator. Moreover, by (1481) . Y*Y = X + 
X'^ < (1 + ||X||)X, and hence h^/‘^Y*Yh^^‘^ is also a trace class operator. 
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Consequently, k*Y = J*kh J is a trace class operator because 

kh~^^‘^ and h^^'^Y* are Hilbert-Sclimidt operators. □ 

Now we are ready to finish 

Proof of Theorem [H From Lemma [9l we know that the quadratic form H 
defines a self-adjoint operator, still denoted by H, such that 

> -^Ti{kh-^k*). 

Let V be as in Theorem [TJ Let 'L be a normalized vector in Q defined in 
([B]). Consider 'L' := Uy'L. From 
that 


7^/ 
a^i 

Moreover, we have 


a 




1 + e/7\I/' J* 


= V 


and cni), it is straightforward to see 

(60) 


7^ 


a:i 


1+J7V1, J* 


V. 


V* 


0 0 
0 1 


V = 


Y 


1 + JXJ* 


(61) 


with X = V*V and Y = JV*JU = JU*J*V. From (IBHIl and (IBTI) . we obtain 


7 ^/ - X a\,, - Y* 
avp/ - y J( 7 My/ - X)J* 


7^' 

a^,! 1 + 


X Y* 

Y 1 + JXJ* 


= V* 


= V* 


7'P 

7^ 


Oiv\ 


1 + J^^J* 


v-v* 


0 0 
0 1 


V 


a 




V. 


Q;\]> J^'^J* 

Recall that 7 >ir and a-q, are finite-rank operators because 'if £ Q. Therefore, 7 >ir/ 
X and a\i,i — Y are also finite-rank operators. Using the cyclicity of the trace 
we find that 

- X)h^/‘^) + mY{k*{a^, -Y)) 


= -Tr 
2 


= - Tr 
2 


= - Tr 
2 


A 

AV 


-X - y* 

a^^-Y J{-f^,-X)J* 




V 


1 


= - Tr 
2 




7-^ 


a:. 


a 




Ot^i J'^^iJ' 

= Tr(^7vt) = ('I/,dr( 0 'L) 


7 ^ a,3 

f, 0 \ / 7^, 

0 JfJ* J Y J^qiJ* 

Thus by the quadratic form expression ([5]), we have 

(T,UvHU{;^') = +ilTv{k*aqr) 

= - X)h^/^) + ii:TT{k*{a^, -Y)) 

+ Tr{h^/^Xh^/^) + 3?TV(ry) 

= {if, dr(e)^') + TT{h^/^Xh^/^) + 5ft Tr(r y) 

for all if £ Q. Recall that we have proved in Lemma fTTl that hfl’^Xhfl'^ and 
k*Y are trace class operators. Hence, 

UvHUt; = dr(^) + Yr{h^l‘^Xh^l‘^) + 5ftTr(ry). 
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Since dr(,f) has a unique ground state |0) with the ground state energy 0, 
we conclude that H has a unique ground state 'ho = Uy|0) with the ground 
state energy 

inf(T(H) = Ti{h^/‘^Xh^/^) +^Ti{k*Y). 

Finally, using (IHHll and (iGTIl we find that = X and = Y. □ 

6. Appendix 

In this appendix we prove two abstract results in functional analysis. 

Proof of Lemma O Let us first consider the case when 9) is finite dimen¬ 
sional. Since = C*C is non-negative, it has an eigenvalue /U > 0 with an 
eigenvector tt / 0. Since {C + y/]I){C — y/]I)u = 0, we have either Cu = y/JIu 
or Cv = y/JIv with v = i{C — y/Ji)u / 0. Thus C has an eigenvector ui E f) 
with a non-negative eigenvalue. Since C = C*, it leaves the orthogonal 
subspace {rti}"'" invariant. By the previous argument, C has an eigenvector 
U 2 € {wi}"*" with a non-negative eigenvalue. By iterating this process, we see 
that C has an orthonormal eigenbasis for i) with non-negative eigenvalues. 

Now we consider the case when Sj is infinite dimensional. Assume that 
is compact. Then by the spectral theorem has an orthonormal eigenbasis 
for t). Note that C leaves every eigenspace of invariant. Therefore, we 
can apply the result on the finite dimensional case to every eigenspace of 
except Ker(C'^). Moreover, Ker(C'^) = Ker(C'*C') = Ker(C). The desired 
conclusion follows immediately. 

Now assume that = 1. Take an arbitrary trace class operator K > 0 
on Sj and consider the operator Ki := K + C*KC. Note that Ki is also 
trace class and Ki > 0. By the spectral theorem, Ki admits an orthonormal 
eigenbasis for Sj and every eigenspace of Ki is finite dimensional. On the 
other hand, CKi = KiC because C = C* and = 1. Therefore, C leaves 
every eigenspace of Ki invariant. Then applying the result on the finite 
dimensional case for every eigenspace of Ki we get the desired conclusion. 

□ 

Proof of Lemma . By taking the expectation of against a vector u 
and optimizing over e > 0, we obtain 

\{u,Lu)\ < {u,Ku)^^‘^\\u\\, Mu. (62) 

Let us show that L is a compact operator. It suffices to show that Lun 0 
strongly for every sequence ^ 0 weakly. We first consider the case 
L > 0. Since ^ 0 weakly, ||nn|| is bounded by the principle of uniform 
boundedness. Moreover, y/Kun —>■ 0 strongly because y/K is compact. 
Therefore, from (|62p we get 

\\LUn\\^ < ||VL|P • \\\^Un\\^ < \\y/L\? ■ \\\^Un\\ ' ||rin|| 0. 

In the general case, when L is not necessarily positive, from (1451) we have 

2L1(L > 0) < e"^l(L > 0)A:1(L > 0) -k e, Ve > 0. 

Since 1(L > 0)iLl(L > 0) is trace class and L1(L > 0) > 0, we conclude 
that L1(L > 0) is compact. Similarly, L1(L < 0) is compact, and hence 
L = L1(L > 0) -|- L1(L < 0) is compact. 
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Finally, since L = L* is a compact operator, it has an orthonormal eigen- 
basis {vn} with real eigenvalues. Using (1621) again we obtain 

Tr(L2) = Y,\{Vn,LVn)\^ < Y,^Vn,KVn) = TV(i^). 

n n 

□ 
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